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We calculate the A polarization in an exact analytical, rotating model based on parameters 
extracted from a high resolution (3+l)D Particle-in-Cell Relativistic hydrodynamics calculation. 

The polarization is attributed to effects from thermal vorticity and for the first time the effects of 
the radial and axial acceleration are also studied separately. 
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I. INTRODUCTION 

In high energy peripheral heavy ion collisions there is a 
substantial amount of initial angular momentum directly 
after the Lorentz contracted nuclei penetrate each other. 
The formed Quark Gluon Plasma locally equilibrates, the 
shear flow leads to local rotation, i.e. vorticity, and then 
it expands, while its rotation slows down. 

Due to the finite impact parameter, the initial stages 
(IS) have a non-vanishing angular momentum [Ull]. For 
the initial stages, effective models as the color glass con¬ 
densate (GGG) or Glauber model are used. In general, 
we use experimental data, to construct a possible IS, at a 
given impact parameter for the participant nucleons, and 
their eccentricity. Early studies neglected effects arising 
from the non-vanishing angular momentum, but interest 
increased recently [31]. 

After many decades of refinements hydro dy¬ 

namical modelling became the best to describe the mid¬ 
dle stages of heavy ion collisions at relativistic energies. 
Thus, rotation and its consequences in peripheral colli¬ 
sions were also studied in fluid dynamical models [3iin]- 

We look at polarization in effects arising from thermal 
vorticity in the exact rotating and expanding model El, 
where we are modeling an appropriate time-period of the 
collision [T^] . Special attention has been given to the col¬ 
lective motion, and to extract it from observables which 
could confirm that such descriptions are indeed plausible. 

We calibrate an exact rotating model based on a 
(3-|-l)D fluid dynamical model, the relativistic particle- 
in-cell method (PICK), to fine-tune the initial parameters 
of the rotating and expanding fireball |12) . 

In Ref. [T3] the differential Hanbury Brown and Twiss 
(IIBT)-method was used to detect rotation in heavy ion 
collisions. 

Without at least some viscosity and/or interaction one 
could not generate rotation from the original shear flow. 
On the other hand to develop instabilities or turbulence 
the viscosity should be small, so that the ratio of shear 
viscosity to entropy density, rj/s, should be of the order 
of h/AirkB, which can be achieved at the phase transition 
between hadronic matter and QGP |14) . 

Thermal vorticity arises from the flow velocity field 
[H], and the inverse temperature field present in heavy 


ion collisions, and it arises mainly due to a non-vanishing 
angular momentum in the initial stage. 

Fluctuations in the transverse plane can generate sig¬ 
nificant vorticity, but in peripheral collisions the initial 
shear flow leads to an order of magnitude larger vortic¬ 
ity. m- This vorticity may be further enhanced by the 
Kelvin Helmholtz Instability (KHI). 

In our formalism, the dynamics of the system af¬ 
ter local equilibration is computed using the relativistic 
(3-|-l)D fluid dynamical model, PICR. This fluid dynam¬ 
ical (FD) computation with small viscosity shows en¬ 
hanced collective rotation due to an evolving KHI. In 
Ref. |16j a simple analytic model for this phenomenon 
is explored using a few material properties: the surface 
tension between the colliding nuclei, the viscosity, and 
the thickness of the flow layer. This enables a classical 
potential flow approximation, in which one may study 
the dynamics of an onsetting KHI. 

A more recent calculation of the onset and effects of the 
KHI is performed in Ref. [H], in which the calibration 
of the exact model takes place. Here, it is pointed out 
that this feature - the enhancement of rotation - is a 
dominant aspect of the (3-|-l)D fluid dynamical model, 
but it is also seen in UrQMD m 

At high energy collisions, we need an initial state 
model, which describes the dynamics until local equi¬ 
libration is reached. There are several options for de¬ 
scribing this pre-equilibrium dynamics, using Color Glass 
Condensate (CGC) helds, parton (or hadron) kinetic the¬ 
ory, or one dimensional Yang-Mills field (or flux tube) 
models [Ull]- In the (3-|-l)D PICR fluid dynamical model 
that we use as our guidance for the FD development, this 
last choice is used. 

It is important to mention that for peripheral collisions 
the initial shear and sometimes even the angular momen¬ 
tum are neglected, while realistic initial state models in¬ 
clude these features [313 HD- 

^From the initial shear flow, in the (3-|-l)D PICR fluid 
dynamical model the general rotation develops gradually 
in 1 to 2 fm/c time. Thus, the Exact model is applica¬ 
ble from this point of time on [T^]. At the energies we 
discuss, by this time the matter is in locally equilibrated 
QGP phase, and the local vorticity develops also. Due to 
the spin orbit interaction the local vorticity and the spin 
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of quarks equilibrate. The essential part of the dynami¬ 
cal development of flow (and other collective mechanical 
processes) takes place in the QGP phase, which is indi¬ 
cated by the Constituent Quark Number scaling of the 
flow harmonics. 

This most significant middle stage of the reaction can 
be modeled by the ’’Exact” model m- The model is 
based on a set of scaling variables, 

( 1 ) 


in terms of the transverse and axial coordinates, x z, and 
y] and the characteristic radius R and axial length, Y, 
parameters. The scaling parameter s = + Sy is also 

introduced, being the scaling variable as it appears in the 
thermodynamical relations. Here we have interchanged 
the y- and z-axes to resonate with choice of axes in heavy 
ion collision literature, in which the reaction plane, in 
which the system rotates, is spanned by and e^, leav¬ 
ing the axis of rotation to be defined by By. 

Ref. [12] calibrates the parameters of the Exact model 
to the (3-|-l)D fluid dynamical model. The parameters 
are extracted for experiments at ^/Smn = 2.76H • TeV 
with impact parameter b = O.lbMax (See Fig. 1 and 2). 
In the (3-|-l)D PICK model, rotation may increase due to 
Kelvin-Helmholtz instability, whereas in the Exact model 
- and the later stages in the experiments themselves - ro¬ 
tation slows due to a transfer of energy to the explosively 
increasing radial expansion of the system. The Exact 
model, therefore, is suited to describe the period from 
the equilibration of rotation up to the freeze-out. 

In El the solution for a flow of conserved num¬ 
ber density, together with a constant, temperature- 
independent compressibility, and a velocity held is de¬ 
scribed. Hence the solutions take form, in cylindrical 
coordinates (r, y, (p ), where r = '/x^~+~z^ with an equa¬ 
tion of motion, r{t) = v{r,t). The Exact model assumes 
a linear velocity prohle both in the radial, r, and in the 
axial, y, directions. This leads to a how development 
where a huid element starting from a point (ro,yo, po), 
and at a later time, t, reaches the point 


r{t) = To 


y{t) = yo 


m 

yjt) 

YitoY 


p{t) = po + 


dtuj{t), 


( 2 ) 


showing explicitly how the solutions evolve in time, ro¬ 
tating and expanding huid. These equations follow the 
time-evolution of the scaling variables in the radial and 
axial directions. This is a Cylindrically symmetric setup 
with X{t) = Z{t), + Z'^it) = R{t) and, in gen¬ 

eral, Y{t) ^ R{t). 

We have chosen the x, z-plane as our plane of rotation, 
with y being the axis of rotation. Our initial angular mo¬ 
mentum, then, points in the negative y-direction, with 


an absolute value of approximately 1.45 • In an at¬ 
tempt to determine new observables, we propose a search 
for A polarization. Although, the polarization could be 
described similarly for all Fermions, we chose the As, be¬ 
cause it is straightforward to determine its polarization 
from its decay to p and tt, (where the p is emitted into the 
direction of the polarization). Actually such an experi¬ 
ment is already performed at RHIC, but the results were 
averaged for A-emissions to all azimuths, while we pre¬ 
dict significant polarization for particles emitted in the 
±x direction in the reaction plane |18j . 

Our expectation is that this polarization, at least in 
part, will be able to account for the polarization as ob¬ 
served in peripheral regions in the first 10-15 fm/c fol¬ 
lowing the impact in a heavy ion collision. 

In order to evaluate the polarization in the Exact 
model we use the parametrization of the Exact model 
based on the realistic (3-|-l)D PICR fluid dynamical cal¬ 
culation [12] , and use the vorticity calculated in the Exact 
model with these parameters in Ref. m- 


II. FREEZE-OUT AND POLARIZATION 

Polarization of As was subject to theoretical studies 
before, both in p+p and in heavy ion reactions. In single 
p + p collisions forward production in small-transverse- 
momentum fragmentation was theoretically studied and 
also observed. These reactions did result in much higher 
polarizations up to about 30% [20] . 

To apply this approach to heavy ion collisions is a com¬ 
plex theoretical problem because several microscopic pro¬ 
cesses can contribute to polarization and these can be 
combined with different hadron formation mechanisms 
pn 122]. In Ref. [21] it was contemplated that the fi¬ 
nal heavy ion results are dependent on the hadronization 
mechanisms, and the effect of the decay products of the 
polarized hyperons on the V 2 flow harmonics, V 2 , were 
studied. Ref. [22] has also studied the sensitivity of A 
production on the coalescence or recombination mecha¬ 
nisms of the hadron formation. 

As the previous works discussed a wide variety and 
complexity of the microscopic description of hadroniza¬ 
tion and the resulting polarization, we have followed a 
simpler statistical picture, based on some simple assump¬ 
tions of a dilute gas of particles, on the ’’Relativistic dis¬ 
tribution function for particles with spin at local thermo¬ 
dynamical equilibrium” [23] . 

This work does not address the mechanisms of 
hadronization and the change of polarization during this 
process. It also barely discusses the equilibrium between 
particle polarization and local rotation in thermal equi¬ 
librium for dilute gases. Thus, this approach is primarily 
applicable to the final hadronic matter. 

We follow the same reaction mechanism as used in all 
(3-|-l)D PICR publications since 2001. We do not as¬ 
sume a 3 stage fluid dynamical process in QCP phase, 
mixed phase and hadronic phase because the fastest adi- 
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abatic development in the mixed phase would take 30-50 
fm/c |24| . Such a long expansion time would contradict 
all two particle correlation measurements showing a size 
and timespan at FO of less than 10 fm. Furthermore it 
would also contradict to the observed Constituent Quark 
Number Scaling and to the observed large Cl abundance. 
The only way out of these problems is supercooling in 
the QGP phase, followed by rapid hadronization PSIES] . 
and almost immediate freeze-out. 

Thus in the PICK fluid dynamical calculations we dis¬ 
cuss exclusively the QGP phase, even for supercooled 
QGP. Based on the mechanical equilibrium, evidenced 
by the Constituent quark number scaling, we have rea¬ 
son to assume that during the FD evolution there is am¬ 
ple time to equipartition the local rotation among all 
degrees of freedom in QGP, due to the spin-orbit interac¬ 
tion. As this is a strongly interacting form of matter the 
kinetic approximation as a dilute gas is not necessarily 
applicable, and the energy momentum and local angular 
momentum should also be carried by the fields. [J We 
have to assume that the rapid hadronization maintains 
equipartition among all degrees of freedom carrying an¬ 
gular momentum. So, based on this assumption we use 
the approach of [23]. 

Actually the same applies the statistical and thermal 
equilibrium among (most of) the abundances of final 
hadron species. This can be understood based on the 
fact that the statistical factors are the same in rapid for¬ 
mation of hadrons as in thermal equilibrium. 

We use the same assumptions for the Exact fluid dy¬ 
namical model as we used for the (3-|-l)D PICK fluid dy¬ 
namics. Based on the above, in the Exact model the en¬ 
ergy weighted thermal vorticity was calculated [Hj. We 
explored the total energy of the system and the energy of 
expansion, rotation, and internal energy components and 
their time dependence. We observed the transfer of en¬ 
ergy from rotation to expansion, hence the rotation slows 
as the system expands until the freeze-out. 

According to the quantum-field-theoretical approach 
[23j , the expectation value of A polarization in an inverse 
temperature field, /3^{x) = u^{x)/T{x), is 

(n^(a;,p)) = ^e;xp<Tr(l - nF)5^/3‘"(a;) —, (3) 

o ui 

where e^po-r is the completely antisymmetric Levi-Civita 
symbol, np is the Eermi-Jiittner distribution for spin-1/2 
particles ((1 — n^) is the Pauli blocking factor), and p is 
the A four-momentum. We integrate this over some vol¬ 
ume, and ultimately over all of space, weighted by the 
number density, normalized by number of particles in 
that volume, leaving a momentum-dependent polariza¬ 
tion four-vector in the participant frame of reference 


^ If we would consider only three valence quarks in kinetic equi¬ 
librium according to [23], then the polarization of a coalesced 
Baryon would be Ilg ~ (llg)^, which would not be measurable. 


^ f J dExp^nFix,p){l - nFix,p))dPf3'^ 

n^[p) - / dEAP^n^(a;,p) 

(4) 

Note that, as opposed to electromagnetic phenomena, 
in which particle and anti-particle will have anti-aligned 
polarization vectors, here it is shown that A and A po¬ 
larizations are aligned in vorticious thermal flow fields. 

While the average values of polarization may be as low 
as 1-2%, consistent with RHIG bounds, in some regions 
of momentum space we see a larger polarization, about 
5% for momenta in the transverse plane and up to a mo¬ 
mentum of 3 GeV/c. Kelvin-Helmholtz instabilities may 
further enhance rotation, hence the thermal vorticity, de¬ 
fined as 

(5) 

and thereby the signal strength increases by 10-20%. At 
LHC energies, there may be 5% A polarization due to the 
corona effect, single nucleon-nucleon collisions occurring 
outside of the reaction zone of the collision itself. So 
attempts should be made to further the understanding 
of this background, and remove it from measurements 
in order to further isolate the A polarization as it arises 
from the collision itself. 

The A polarization is determined by measuring the an¬ 
gular distribution of the decay protons in the A’s rest 
frame. In this frame the A polarization is no(p), which 
can be obtained by Lorentz boosting the polarization 
n(p) from the participant frame to the A’s rest frame, 

ini, 

n„w = n(p)--5|J^n(ri.p, (6) 

where (p*^,p) is the A’s four-momentum and m its mass. 

Based on this equation we see that in order to maxi¬ 
mize polarization, we need to choose momenta for the A 
such that they lie in the reaction plane, hence we fix p 
in the positive x-direction. 

III. SOLUTION FOR THE A POLARIZATION 

As the A is transversely polarized, H^Pp = 0, one can 
confine himself to the spatial part of H^. The simplified 
spatial part of polarization vector is: 

^ he J dVnF{x,p) (V x /3) 

^ 8m J dVnpiXjp) 

J dVnF(x,p) (dt/3 + V/3°) 

8m J dVnF{x,p) 

where npix^p) is the phase space distribution ofthe As. 
In a previous calculation m. the p dependence of riF, 
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FIG. 1. (color online) The direction of axes, as well as the 
momentum, p, and flow, /3, vectors. The azimuth angle is 
measured from the direction of the p-vector, i.e. from the 
i-axis. 


was considered negligible in the integral and the time 
derivative and gradient terms were also assumed to be 
smaller. The present calculation shows that in general 
these terms are not negligible and that which terms are 
dominant depends on the particular conditions. 

We adopt the parametrization of the model from Ref. 
m. with the initial conditions Rq = 2.5 fm, Yq = 4.0 fm, 
Rq = 0.20 c, Yq = 0.25 c, ojQ = 0.1 c/fm, k = 3/2, Tq = 
300 MeV. For this configuration Etot = 576 MeV/nucl. 


A. The denominator 


We first perform the integral in the denominator: 

R +Y 271 - 

A{p) = j dV Up = J r dr J J d(j) np^Xjp) . (8) 


-F 0 


According to Eq. (3) in Ref. [19] in terms of the scaling 
variable, s, we have: 


n = no — v(s) , 

, , 1 /I 

iy(s) = —- exp 
r(s) 


du 

t{u) 


= 1 • exp — - / du 


(9) 


( 10 ) 


where the simplifying choice of r(s) = 1 is used in the 
last step. Therefore: 


/ ^ bb -is 
n{s) = no—e = . 


( 11 ) 


The EoS is assumed to be: e(s) = KT{t)n{s) and the 
energy density e(s) is calculated as in Eq. (29) in Ref. 


m, therefore: 


n(s) = 


KT{t) 


Cn 

kT 


e 


2 e 


Sp 

2 


( 12 ) 


where Cn = ^■ 

^From Ref. [IH], the Fermi-Jiittner distribution is: 


nF{x,p) 


1 


1 

gP'^/3,7-5 


eCT 


(13) 


where the ^ = /i/T, and p, is the chemical potential. The 
thermal flow velocity, = u^(x)/T, is different at 

different space-time points x. 

The invariant scalar density for the Jiinttner distribu¬ 
tion is: 


n - 

{2phY "Co 

where the C/"^ = A'Km?TK 2 {m/T)/{2TTfiY. With Co and 
n{s) = n, the Fermi-Jiittner distribution can be written 
as: 


eP/T 


Con(s) 


(15) 


Now we introduce cylindrical coordinates for the location 
in the configuration place x = (r,y,(j)), and using the 
scaling expansion model [nma with the scaling vari¬ 
ables s , Sr ,Sy. Now, substituting Eqs. PPl ) into the 
denominator of n(p), and parametrizing the range of in¬ 
tegrations as in [ini one obtains: 


A{p) = 


CnCq 

kT 


naY 

dy exp ( - 


l-aY 

p27r 


yY\ 
2^2 ) 


d(j>e 




l<bR 

rdr exp 


'0 


2R2 


(16) 


The scalar product in cylindrical coordinates takes the 
formp'"/?^ = (p°,p)(/3o,/3) = p°Yo -p(3 = P°Yo -PrPr - 

Py/dy ~ P4>fdip- 

In our integral the p^ is given or ’fixed’ as the argument 
of n(p), while the /3 = f3{x) is changing. The integra¬ 
tion with respect to (j) starts from the direction of the 
p-vector. According to the Eq. (5) in [T^ : 

V = VrCr + VyBy = ^TBr + OJTB^ + yUBy, and 

/3 = u'^/T — ^vjT. Thus in the integral for we exploit 
the fact that in the Exact model the radial, r, and axial, 
y, components of the thermal velocity, /3, do not depend 
on (/, while the tangential component does not depend 
on y, i.e. = jrix/T, but its direction is changing 
with respect to the direction of p. As the integral is over 
the whole 27r angle we can start it at any point of (j), so 
we start it from the externally given p-direction. Con¬ 
sequently, with this choice of the x-axis, p = (jpr,Py,0), 
and Pz = Ptp — 0. In this azimuthally symmetric, exact 
model it is sufficient to calculate n(p) for one direction 
of p in the [x, z]-plane. 
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The direction of the thermal flow velocity, /3, is tan¬ 
gential to the direction (p, i.e. it points to the e 0 +,n./ 2 - 
direction. Thus the scalar product is: 

p-(3{r, y,(t>) = \px\l3r cos {4>)+pyl3y + \p^\l3^cos[(j)+'^ , 


where (p is the azimuth angle of the position around the 
y, rotation axis, counted starting from the x-axis. See 

Fig. [2 

So, inserting the last expression for p^Py into the last 
term of the integral Eq. (16), the integral with respect 
to (p will take the form: 

J = J = 27r/o(\/a2-t&2) , 

(17) 


where a = \px\Pr = \Px\ 7 Rr/TR and b = \px\P 4 , = 
\px\ 7 ruj/T, and we used integral no. 3.338(4) in [57]. 
If we define 


C 3 = 


/ PxiR \ 

, (Pxl^^P 

[ TR ] 

+ V T ) 


^ ^ \Px\l 

T 


\I{R/Ry 




then y/a? + h'^ = csr, and: 


p27T 

/ d(pe~^‘"^>^ = IT" ^PvPy X 2ttIq {c^r). (18) 

Jo 

Now, substituting this back into Eq. (16): 

Mp) ^ j nF{p,s) = 


CnCq 

kT 


aY bR 


I dr, J rdr ^ 


-aY 0 


X ■ 27rJo (cg r) . 


(19) 


Now we may use the same simplifying non-relativistic 
assumption as in Eq. (5) of Ref. |19|. i.e. we approximate 

by as V = VrBr + VyBy + v,pe,f, = ^rBr + yyey + 
urre^, and thus 7 = 1 . It follows then: 


A{p) 


CnCq 

kT 


dV npip, s) 


aY 


27re J ex.p(ciy-C2y‘^)dy 

-aY 

bR 

X J r/o(c3r) exp(-C4r^)dr , (20) 


where Ci = PyYj{YT), C 2 = l/(2y^), C 4 = l/(2i?^) are 
constants. 

Now we assume an infinite system with scaling Gaus¬ 
sian density profile, so that the integrals are evaluated 



pJGeV/c) 


FIG. 2. (color online) The polarization of A particles, ni(p), 
in the participant Center of Mass (CM) frame for the first 
term containing the (V x/3)-contribution, at time t — 0.5 fm/c 
after the equilibration of the rotation, in the Exact model. 
The polarization, IIi(p), points into the —y-direction and 
changes from —1.5% at the CM-momentum {px = Py = 0), 
to —8% in the corners, in 1% steps per contour line. The 
negative percentage indicates that the polarization is in the 
—y-direction. The structure is just like that of the energy 
weighted vorticity. Due to azimuthal symmetry of the Exact 
Model the px and pz dependence of II are the same. 


up to infinity, i.e. the parameters a = 00 , b = 00 . Thus, 
the y component integration in Eq. (20) is calculated as: 
/*+00 


^ciy-C 2 y 


dy = . — exp 


Led' 


( 21 ) 


where we used the integral formula No. 2.33(1) in 
and erf (-I- 00 ) = 1, erf (— 00 ) = —1. 

For the integration of r component: 

r+°° 2 

/ r/o(c 3 r)e dr 


1 


I'A'i 


I ^3 


C3^/bI I 8 c4 ) 

where the M_y ,z{z) is the so called ’Whittaker Function’, 
No.6.643(2) in ’nil . 

Now, we obtain the final form of Eq. (21): 


A{p) = 


_ 27rv^ CjvCo yO/T ( _Ci_ 


RT C^y/OzCA 


X exp 


f —1 

V 8 C 4 / 


M 


-i 0 


exp 

(—)■ 

V4c4/ 


( 4 ) 


(23) 


However, in the relativistic case, the integrations 
with respect to y and r can not be performed ana¬ 
lytically, because of the presence of the factor 7 = 

\lJx-vl-vl-vl . 
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B. The numerator 

Ref. [19] calculates the energy weighted vorticity, 
which is azimuthally symmetric, i.e. independent of the 
azimuthal angle 4 >- In. ths dehnition of the polarization, 
Eq. §, we have npiPyx) = enpiPjx) for As with 
momentum p. In [19j . however, the energy weighting 
is performed with the total energy density of the fluid 
Etot = Eint + Ekim which in general is not the same as 
enpiPjx). On the other hand the bare vorticity is just 
a constant in the non-relativistic Exact model, while the 
EoS may be more general and it may lead to more in¬ 
volved R{t) and Y (t) dependence than the ideal Jiittner 
gas approximation would allow. 


mentioned above (beginning of sec. Ill), the thermal 


vorticity is h(V x /3) = —0.13 at t = 0.5 fm/c, and it 
decreases very slowly with time, about 1-2% per 1 fm/c. 
This constant vorticity will make the numerator simple: 


B{p) = J dV nF(V X P) = X A{p) 


(25) 


Therefore, the first term of polarization vector, i.e. Eq. 
Q will be: 


ni(p) 


he J dVnF(x,p) (Vx/3) ^ heu 
8m J dVnp(x,p) 



FIG. 3. (color online) The absolute value of A-polarization, 
n2(p), in the participant Center of Mass (CM) frame for 
the second term containing the (5t/3)-contribution, at time 
t = 0.5 fm/c after the equilibration of the rotation, in the 
Exact model. The polarization changes from zero at the 
CM-momentum (p^, = py = 0), up to 20% in the corners 
at px = —4GeV/c, in 2.5% steps per contour line. In the 
corners at = 4GeV/c, the polarization is 12%. This 
second term is orthogonal to p, and it is smaller, especially 
at CM-momenta, where it is negligible. This term arises 
from the expansion, which is increasing rapidly in the Exact 
model with time and also increases with the radius. At large 
radius the larger expansion leads to larger momenta. The 
structure of the 2nd component of polarization arises from 
the asymmetries of the different components of 112 (p) 


which means the polarization vector arising from the vor¬ 
ticity, ni(p), in the Exact rotation model is a constant, 
(although time dependent), and parallel to the j/-axis. 



P,(GeV/C) 


FIG. 4. (color online) The x component of the A- 
polarization, n 2 a:(p), in the participant Center of Mass (CM) 
frame for the second term containing the (dt/3)-contribution, 
at time t = 0.5 fm/c after the equilibration of the rotation, 
in the Exact model. The polarization vanishes at the CM- 
momentum {px — py = 0), and changes from zero up/down 
to ±8% in the corners, in 1% steps per contour line. This 
term arises from the expansion, which is increasing rapidly in 
the Exact model with time and also increases with the radius. 
At large radius the larger expansion leads to larger momenta. 


Thus we use the direct, non-relativistic vorticity val¬ 
ues, U!(t), from Ref. [121, and not the presented energy 
weighted vorticity. I.e. 

V X 13 = -2u}{t)ey/T{t) , (24) 

so that the thermal vorticity has only y-directed compo¬ 
nent in the Exact model. With the model parameters 


One may add the Freeze-Out (FO) probability to the 
integral. According to the Ref. [2H[, the FO probability 
is Ws = (Pfi^s) (P ■ '*^(®))i where the approximation is 
used that the FO direction, is parallel to the flow ve¬ 
locity, u(x) = "fv^x). In the first term of the numerator, 
which depends on the constant y-directed vorticity this 
FO probability influences the numerator and denomina¬ 
tor the same way, so the effect of the two integrals cancel 
each other in the FO probability also. 
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C. The second term 


The numerator in second term of polarization vector 
reads: 

C{p) = I dVnp{x,p) {dt(3 + VI3°) . (27) 


If, in the non-relativistic limit, 7 = 1 is assumed, then 
V/3° = 0, and dtP = dt{v/T), so we have to evaluate only 
the first term of the sum in the integrand. According to 
Ref. [El [29], the time derivatives of velocity are: 


OtVr 






■R 



' W 

ch + 2 ^u; 

'Y 

y 2 - 

Y ~ 




y = C7y, 
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(28) 


where C 5 = (i?/i? — Ft?/— w^), cq = {lu + 2{R/R)uj), 
and C 7 = {Y/Y — Y"^jY"^). 
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FIG. 5. (color online) The y-component of A-polarization, 
n2(p), in the participant Center of Mass (CM) frame for 
the first term containing the (dt/3)-contribution, at time t = 
0.5 fm/c after the equilibration of the rotation, in the Exact 
model. The polarization changes from zero in the middle to 
—8% at px = ±4GeV/c, in 1% steps per contour line. This 
y-component points into the axis-direction just as the first 
term. Hi, thus these two are additive. The ^-component of 
112 (p) does not depend of Py , as shown in Eq. (361. 


Therefore, 


dt(3 = (c 5 rer + cere^ + cpyey) jT , 

provides the time-components of the vorticity in the 
three spatial directions. Here, as the model is symmet¬ 


ric, dtPy vanishes, and with the model parameters men¬ 
tioned above (sec. Ill), at t = 0.5 fm/c and r = I fm 
^dtPr = 0.024 and = 0.009. Both these vorticity 

components decrease slowly with time by about 0.0005 
in I fm/c. 

Eq. (27) is a volume integral of a vectorial quan¬ 


tity, which is not convenient to perform in cylindrical 
coordinates. So we transform it into Cartesian coordi¬ 
nates: By = cos (f>ex + sin = — sin (pe^Y cos p . 

Therefore, T ■ dtf3 = (05 cos p — c% sin p')rexY (05 sin p + 
C6 cos p)re^ + cry By 


The integral of Eq. (27) can be expanded as: 


C'(P) = J dVnF{x,p)dtP 


_ CjyC'o po/T 


kT 


rdrdpdy exp (ciy — C 2 y^) 


X exp (^acos p — bsinp — c^r'^) dtf3 , (29) 

where a and b are defined after Eq. (nl 
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FIG. 6. (color online) The y-component of A-polarization, 
II(p), in the participant Center of Mass (CM) frame for the 
second term containing the (9t/3)-contribution, at time t = 
0.5 fm/c after the equilibration of the rotation, in the Exact 
model. The polarization is -1.5% at the CM-momentum (pa, = 
Py = 0), it is —16% in the corners. The change is in steps of 
2% per contour line. 


It is convenient to dehne an integrating operator, A as: 

A = J dVnpiXjp) X 


r dr dp dy e 


Ciy—C2y a cos (p—b cos (f>—C4r'^ 


and then Eq. (29) will be: 


C{p) = A ■ dtf3 = ^ {Ib^ + Jb^ + HBy) , (30) 











































where we defined: 


I = A ■ (c 5 cos 4> — cq sin (ji) r , 

J = A ■ {c^ sin (j) + Cq cos 4>) r , 

H = A ■ C’jy . 

Using the integral formula No. 2.33(6) of [27] the func¬ 
tion H becomes: 

^ _ 27rv^C'jvC'o ^_p^/j C 7 C 1 

kT 2C3C2-yC4C2 

■ (31) 

The function I can be expanded as a function of in¬ 
tegrals over (/), r and y. The integral over cj) brings in 
the Bessel function, 2 TrcsIi{c 3 r)/C 3 (See No. 3.937 (1) 
and ( 2 ) of HZ]), where Cs = (csa' — cgb'), a' = a/r = 
\px\R/TR, and b' = b/r = \px\u;/T. Subsequently, the 
integral with respect to r brings in the ’Whittaker Func¬ 
tion’ and then the final form of I after performing the 
separable integration with respect to y leads to: 

j _‘^'!^\ArCNCo cg 

kT c|c 4 y(c^ 

x“KS)“H4|)"-M(ii) ^ (32) 
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FIG. 7. (color online) The a-component of A-polarization, 
n2(p), in the participant Center of Mass (CM) frame for 
the second term containing the (dt/3)-contribution, at time 
t = 0.5 fm/c after the equilibration of the rotation, in the Ex¬ 
act model. The polarization vanishes at the CM-momentum 
{Px = Py = 0), it is ±3% in the corners. The change is in 
steps of 0.5% per contour line. The corners at py = —4GeV/c 
are positive while at py = 4GeV/c are negative. 


Evaluating the integral J is similar to J: 


J = 


^27ry(7rCArC'o -p^/T 


cg 


kT 


C 3 C 4 A/C 2 


xexp(i)exp(i)M-.,.(i) , (33) 


,2 

4c4 3 


where the only difference is: Cg = {c^b' + c^a') compared 
to Cg in /. 

Then, substituting /, J, H back into Eq. (30), one can 
obtain the analytical solution for numerator in second 
term of polarization vector as: 

C(p) = J dVnF{x,p) dtl3 = ^ {lex+Je^+Hey) 


kT^ 

C 8 


r .2 

V 8 C 42 ' 


exp 

Cg 




C3C41/C2 ’ 2 MC 4 '^ C3C41/C2 ’2V4C4/ 


- 

2 C 3 C 2 VC 4 C 2 


(—)' 

V4c4/ 


Dividing this by A{p), i.e. Eq. (23), one gets: 
C 8 


Cg ^-l.i 


C{p) 1 

A{p) TlcsPci M_i^ ' C3y(cIM_i o' 


(34) 


C 7 C 1 
2 c 2 

( 35 ) 


Then, we obtain the second term of polarization vector: 


n 2 (p) 


ftp ^ Cjp) 
8 m A{p) 


ft 

PyCg 

|Pa:|cg ^-l.i ^ 

8mT 

C3^M_1o 



, f\Px\c7Ci 

PyCs 


+ V 2c2 



(36) 


As we can see, and as is given also by the definition, 
Eq. (0, the second term of polarization is orthogonal to 
the particle momentum: 


n2(p) T p , (37) 

thus if we use the choice that p should be in the [x,y]- 
plane and its z-component should vanish, then the y- 
component of n 2 (p), should depend on px only, see Fig. 

[51 


IV. THE FREEZE-OUT STAGE 

The fluid dynamical model is in principle not adequate 
to describe the final, post freeze-out (FO) particle distri¬ 
butions, the abundance of the particle species and also 
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their polarization. This is so because the post freeze-out 
distributions must not be in local thermal equilibrium 
and must not have interactions among the final emitted 
particles. Furthermore, the emitted particles should not 
move back into the interacting zone, i.e. towards the 
pre-FO side of the FO hypersurface. How to handle the 
freeze-out is described in great detail in m- It indi¬ 
cates two ways to handle this process: (i) consider the 
post-FO matter as if it has an Equation of State (EoS). 
This is only possible if the post FO EoS is that of a non¬ 
interacting ideal gas and the FO hypersurface is timelike, 
(ii) The other approach is that the post FO matter is de¬ 
scribed by a dynamical model with weak and rapidly de¬ 
creasing interaction, like UrQMD or PACIAE, matched 
to the QGP fluid on the FO hypersurface. The change at 
crossing this hypersurface is in general significant, as the 
pre-FO matter is strongly interacting, supercooled QGP, 
while the post-FO matter is weakly interacting and has 
different (usually less) degrees of freedom in both situa¬ 
tions. The FO across the hypersurface is stronger if the 
latent heat of the transition is larger. 

The precise way to perform this transition is described 
in [30] . This method is demonstrated in several earlier 
fluid dynamical model calculations (also using the PICK 
method), for precision calculations of flow harmonics. 

As mentioned in the introduction, at high energies 
(RHIG and LHG) the Constituent Quark Number Scal¬ 
ing and the large strangeness abundance clearly indicate 
a supercooling and rapid hadronization. Furthermore at 
these energies the transition is in the crossover domain 
of the EoS, thus the expected changes are smaller, and 
the major part of the EO hypersurface is time-like, which 
allows to use ideal gas post FO distributions, as we do it 
here using the method of |23]. These are the conditions 
which make the changes in mechanical parameters (e.g. 
v) small at freeze-out while the temperature changes are 
larger [50] . 

Thus, just in the case of Constituent Quark Number 
Scaling, we assume that other mechanical processes like 
mechanical polarization will not significantly change at 
freeze-out at RHIG and LHC energies. This conclusion 
is restricted to local thermal and flow equilibrium, and 
should not apply to some of the microscopic processes, 
which dominate p + p reactions. 

Also, in case of freeze-out through space-like FO hyper¬ 
surfaces, the mechanical parameters change significantly, 
the post-FO distribution is far from a thermal distribu¬ 
tion (it is a Cut-Juttner of Canceling Juttner distribu¬ 
tion), and thus the conditions of [23] that we use, are not 
satisfied. 

In this connection we may mention that in earlier re¬ 
lated publications, previous experimental A polarization 
measurements, which were negative, were discussed. It 
was pointed out that polarization as measured was av¬ 
eraged for all A particle directions. Here, as well as in 
the previous detailed PICR fluid dynamical calculations, 
it was emphasized that polarization should be measured 
after finding Event by Event the reaction plane and the 


center of mass of the system. Significant polarization can 
only expected for particles emitted in selected directions. 

Preliminary experimental polarization studies in the 
RHIG Beam Energy Scan program along these lines are 
promising m, and may lead soon to positive quantita¬ 
tive results. At this point of time the present relatively 
simpler EO treatment of the model calculations with con¬ 
stant time FO are sufficient, and can be refined when 
quantitative experimental data will be available. 


A. Conclusion 


Finally, adding Eqs. (361 and (26) we get the analytical 


solution for A-polarization in the Exact model: 


n(p) = 


8mT 
M 


PyCg ^-i.h 


M_i n/ 

2 


C3v^ M_1q 

( 


( 2ew - ^ X 


By + 


\Px\c7C1 

2 c 2 


C3v^ 

PyCs 
C3y/^ M_1 q 


ri), 

in/ 


(38) 


Notice that Eq. 
the non-relativistic limit. 


(38) is the analytical solution in 
The ’Whittaker Function’, 


is the confluent hypergeometric function. For 
the relativistic case, the integrations of the A-polarization 
vector cannot be performed analytically, because of the 
presence of 7 = 1 /I — v'^ — Vy — which will make 

the integrations more involved. Thus, a numerical solu¬ 
tion for the A-polarization would be needed. 

The effect of vorticity is shown in Fig. The non- 
relativistic Exact model can handle reactions with mod¬ 
est energy and modest rotation, so the overall vorticity 
and the resulting polarization is not too large. Further¬ 
more the rotation and vorticity decrease with time while 
the radial and axial expansion increases. This expan¬ 
sion leads to the second term of polarization, 112 , which 
depends on dt(3 (while the V/3° terms vanishes in the 
non-relativistic approximation). Due to the simplicity of 
the Exact model, the vorticity arising from the shear flow 
of the peripheral initial state is constant in space and de¬ 
pends on the time only. However, due to the construction 
of thermal vorticity, both the angular momentum and 
the temperature in the denominator decrease with time, 
thus V X ,3 is hardly decreasing with the time, and it has 
a significant value, -0.13, in natural units. At the same 
time in this model the time-dependent vorticity is smaller 
by almost an order of magnitude. The time-dependent 
vorticity components also decrease faster than the one 
originating from the initial shear flow. 

Nevertheless the second term in the polarization is of 
comparable magnitude to the term arising from local vor¬ 
ticity. See Fig. 

The presented plots are such that Px points into the 
direction of the observed A-particle, while the Py is the 
axis direction. All results should be either symmetric 
or antisymmetric for a ±py change. On the other hand 
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reversing the Px axis must not change the data, as the x- 
axis is chosen to be the direction of the argument of n(p), 
which must be azimuthally symmetric in the [x, yj-plane. 

The polarization arising from the dynamics of the ra¬ 
dial and spherical expansion, 112 , was not discussed be¬ 
fore in the literature, as the dominance of the vorticity 


effect was anticipated and studied up to now. The 112 
plots in Figs. [Tj show the components of the po¬ 

larization arising from the dynamics of the spherical ex¬ 
pansion. The most interesting y-component arises from 
the x-component of the momentum and the ^-component 
of the thermal velocity change (Fig. [^. 




pJGeV/C) P,(GeV/C) 

FIG. 8. (color online) The (a) radial, x, and (b) axial, y, components of A-polarization, IIq(p), in the A’s rest frame. For 
noa;(p) the contours represent changes of 1% from —9.5% in the upper left-hand corner to 9.5% in the upper right-hand corner, 
whereas the contours of noy(p) change in steps of 2% ranging from Iloy = 0 (!) at the CM momentum (p^ = Py = G) to —12% 
for Px = ±4GeV/c at the edges. Both plots are asymmetric due to the Lorentz boost to the A rest frame. 


Now if we study the axis directed components, this 
is given by liy = IIij, -|- 112^. Both these terms have 
a negative maxima of the same magnitude, —8%, at the 
corners, Px^Py = ±4GeV/c, thus these terms add up con¬ 
structively and result in A-particle polarizations reaching 
-16% at high momenta. At small momenta the polariza¬ 
tion is still the same sign but has a reduced value of the 
order of 1.5% arising from the vorticity (Fig. |^. 

In this Exact model the x and z components of the po¬ 
larization arise only from the second term, 112 (p). The 
X component is reaching ±8%, while the z component 
is smaller, it reaches about ±3%. These both are asym¬ 
metric for ±py change, and show an opposite symmetry. 
The x-component is proportional to Py and the dynamics 
of radial expansion. Thus it follows the signature of Py, 
Fig. El The z-component is proportional to Py and the 
dynamics of radial expansion, thus it follows the signa¬ 
ture of Py, Fig. The z-component is proportional to 
PxPy and inversely proportional to PyPx, Fig. These 
two effects compensate each other so the maxima of the 
polarization are smaller and the symmetry is opposite 


to that of the x-component. This term is sensitive to 
the balance between the axial expansion and the radial 
expansion in the model. 

The A polarization is measured via the angular dis¬ 
tribution of the decay protons in the A’s rest frame, as 
shown in Eq. The resulting distribution is shown in 
Fig. § This new study indicates that the dynamics of 
the expansion may lead to non-negligible contribution to 
the observable polarization. The structure of noy(p) is 
similar to the one obtained in Ref. [TB], but here the 
contribution of the ’’second”, dtf3 term is also included, 
which makes the y-directed polarization stronger at high 
Px values, 12%, while it was 9% in Ref. [18], both in 
the negative y-direction. Furthermore, the second term 
changes the structure, of the momentum dependence of 
noy(p), and it becomes ±px asymmetric. 

Recently the vorticity and polarization were also stud¬ 
ied in two fluid dynamical models [32] . The initial states 
that were used from Bozek and Gubser neglected fully 
the initial shear flow in the central domain of the reac¬ 
tion, in contrast to other models where this is present 
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[D1 El 133ISD- This results in negligible thermal vor- 
ticity in the central domain of the collision (Figs. 3, 
13 of Ref. [33])) and consequently a negligible polariza¬ 
tion from the vorticity from the ’’first term” discussed 
here. Thus, the observed vorticity arises from the ’’sec¬ 
ond term”. 

On the other hand there is qualitative agreement be¬ 
tween Figs. 12 of Ref. |33| and this work in the sense 
that only the y-directed (i.e. [x,z] or [x,y]) component of 
the vorticity leads to an overall average net polarization. 
This arises in both models from the initial angular mo¬ 
mentum and points into the —y-direction. In Ref. |32j 
this arises as a consequence of viscous evolution of the 
initial, vorticity-less flow, while in our Exact model it is 
present in the initial state. 

Recent preliminary experimental results reported for 
the first time |31j . significant A and A polarization for 
peripheral collisions at RHIC for beam energies y/sJjN = 
7.7—39 GeV aligned with the axis direction of the angular 
momentum of the participant system. Furthermore the 


A and A polarizations were pointing in the same direction 
confirming our approach. 

In this work we analyzed and compared the two terms 
of polarization, in the Exact model. Including both ro¬ 
tation and expansion, and vorticity arising from both of 
these effects enables us to study the consequences of the 
two terms separately. This study indicates that the as¬ 
sumptions regarding the initial state are influencing the 
predictions on the observed vorticity, while in all cases 
observable polarization is predicted. 
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